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Abstract — Calculations of NIEL are reported for heavy ions 
prominent in the space environment for energies ranging from 
200 MeV per nucleon to 2 GeV per nucleon. 

i. Introduction 

M ethods have been described for calculating nonionizing 
energy loss (NIEL) due to Coulombic interactions for 
many ions and targets at energies ranging from threshold to 
the GeV range [1-31. For protons and helium ions at energies 
above 100 MeV per nucleon, nuclear non-elastic mechanisms 
also contribute substantially to NIEL. Results that include 
these interactions have been reported [1,4,5]. There still exists 
a large gap in our knowledge of NIEL for relativistic heavy 
ions, as shown in Fig. 1 , which delineates what NIEL values 
are known as a function of incident ion atomic number, Z, and 
incident ion energy for a given target material such as silicon. 
The large area in the Figure corresponding to relativistic 
heavy ions is addressed in this work. 
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Fig 1. Delineation of known NIEL values and those values addressed in this 
work. 
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At energies in excess of about 100 MeV per nucleon heavy 
ions can undergo fragmentation interactions where the 
incident ion and target are reduced to isotopes with lower 
charge and mass. Due to the large number of reaction 
products that emerge from these high-energy heavy ion 
reactions, at first glance calculating NIEL appears to be rather 
prohibitive. Consequently the effects of nuclear fragmentation 
interactions on NIEL have not been reported for heavy ions 
with atomic number greater than that of helium. Since heavy 
ions are present in the space environment (for example see 
ref. 6 for a detailed description), the contribution they make 
to displacement damage in detectors, microelectronics and 
solar arrays remains an open question. This question is 
becoming more significant as missions to the Moon and Mars 
are being proposed in NASA’s new space initiative. On the 
other hand ionization effects of these fragmentation 
interactions have been extensively studied because of their 
potentially damaging biological and material effects [7,8]. 

In this paper we first examine the contribution to NIEL 
from heavy ion fragmentation reactions. This includes 
fragmentation of both the incident (primary) ion and die 
target. These results are obtained using experimental data on 
fragmentation cross sections [9,10] combined with the 
relativistic expression for NIEL [1,2]. Adding these results to 
the relativistic Coulomb contribution of the primary ion leads 
to estimates of NIEL for high-energy heavy ions in the space 
environment. 

II. FRAGMENTATION INTERACTIONS 

A. Primary Ion and Target Fragmentation 

The current picture of the fragmentation process is one 
where a very fast-moving ion collides with a target nucleus. 
Both are reduced to smaller nuclei and a number of other 
particles, primarily nucleons. The primary fragments, i.e., 
those coming from the incident ion continue to move without 
appreciable change of direction or velocity except for a 
relatively small energy exchange with the target [10]. The 
target fragments emerge in an excited state and with a 
significantly lower kinetic energy than the primary fragments. 
Although ionization studies have typically only considered 
the contribution of the primary fragments, here we consider 
the effect of tire target fragments as well. 

Since fragmentation of the incident ion has been studied 
the most, it will be considered first. The production of each 
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primary fragment has a partial cross section associated with it, 
which is a measure of the probability of its production. There 
are generally a large number of primary fragments that 
emerge. For example, an incident ion such as argon that has 
40 nucleons can result in as many as 36 isotopes in the 
fragmentation process. The partial cross section associated 
with the production of a particular fragment species can be 
calculated by first determining the partial cross section for an 
incident proton. An example of a basic equation for this is [9] 

a p =<T 0 f(A)f(E)e~ PAA ex V (-R\z-SA + TA 2 \'')Qtf ( 1 ) 

Here CT 0 is a normalizing factor, while / (A) and / (E) 
apply to heavy targets with Z >30. The factor exp (— PAA) 
accounts for the reduction of cross sections as the difference 
of product and target mass increases. The second exponential 
yields the distribution of cross sections for the production of 
different isotopes of an element with atomic number Z. The 
parameter lQ is related to nuclear structure as described in 
[1 1]. The quantity T] depends on the pairing of protons in the 

nucleus (e.g. even-even, even-odd). Finally, E, accounts for 
enhancement of light evaporation products as described in 
f 12]. Tabulated values of these parameters are provided for a 
range of projectile-target combinations [9, and references 
therein]. 

The cross sections of nucleus-nucleus (N,-N,) 
interactions can now be calculated by scaling the proton 
partial cross sections. The scaling relation is 

0-(N, - N 2 ) = (T p ■ S c ■ £ n £ ( 2 ) 

N, is the ion which is fragmented in a collision with the 
target N,. The factors £ n ,£ L ,£ 1( and £ A are corrections 
for neutron deficient products, light products, nucleon 
stripping, and reactions with large change in mass AA , 
respectively. Tabulations are available for the appropriate 
factors in Eq. (2) as a function of the projectile and the target 
as well as the energy [9]. The above equations permit the 
calculation of partial cross sections for most ion-target 
combinations of interest. 

The partial cross sections determine the probability of a 
particular isotope being produced in a fragmentation 
interaction. In addition to that information the average energy 
of the fragment is needed in order to determine NIEL. A very 
standard assumption used in analysis of ionization data is that 
the nucleons in the fragment have the same velocity and 
momentum as the original projectile [7,8,13]. This is the so- 
called straight-ahead approximation. In the examples shown 
here this assumption is applied. 

In addition to the high-energy fragments of the incident 
ion, the contribution of target fragments to NIEL is also 
considered. When the incident particle is a proton or helium 
ion the target fragments generated contribute to NIEL and the 
process is commonly referred to as spallation. Analysis of this 
proton and alpha particle nuclear inelastic contribution to 


NIEL for a variety of target materials [4,5] shows that it 
scales approximately as the reaction cross section. This 
provides motivation for estimating the inelastic NIEL for 
heavier ions by scaling the alpha particle inelastic NIEL by 
the known ratio of reaction cross sections. It turns out that this 
estimate of the inelastic NIEL is sufficient for our purposes 
because it is much smaller than the Coulombic contribution 
for many incident heavy ions and relativistic energies of 
interest. 

III. Results 

Tables 1 and 2 illustrate the contribution of the incident ion 
fragmentation process to NIEL relative to the Coulomb 
contribution. Table 1 shows results for carbon ions incident 
on a carbon target at energies ranging from 250 to 2 1 00 MeV 
per nucleon. Table 2 shows results for 600 MeV per nucleon 
carbon, oxygen, argon and iron incident on carbon targets. 
The middle column in both tables is the Coulomb 
contribution to NIEL for the incident ion calculated using the 
relativistic form for Rutherford scattering [1,2]. The NIEL 
associated with the primary fragments is also calculated using 
this model. For an incident carbon ion, there are 1 1 isotopic 
fragments ranging from 1 1C to 6Li and 6He. The column on 
the right is the contribution to NIEL from the fragmentation 
of the incident ion. These results were obtained by weighting 
the values of NIEL for each fragment by the corresponding 
partial cross section to obtain the average NIEL of the 
primary fragments. It is seen that the contribution of the 
primary fragments to the total NIEL is very small. This is a 
direct consequence of the large difference between 
fragmentation cross sections and the relativistic Coulomb 
cross sections for the incident ion. The total sum of the partial 
fragmentation cross sections ranged from 240 to 260 mb. In 
contrast the primary carbon ions had Coulomb interaction 
cross sections that ranged from 187 to 983 b. Examination of 
cross sections for a number of different ions show a similar 
effect and this is expected to be characteristic of all 
relativistic heavy ion interactions of interest. 

Table 1 

CALCUIATED NIEL FOR HIGH ENERGY CARBON IONS 
ON CARBON TARGETS 


Energy/Nucleon 

(MeV/n) 

Primary Ion NIEL 
(MeV»cm 2 /g) 

Primary Fragment NIEL 
(MeV»cm 2 /g) 

250 

9.39x10-3 

1.58x10- 6 

600 

4.83X10- 3 

1.55x10- 6 

1050 

3.29X10- 3 

1.56x10-6 

2100 

2.06x10-3 

1.70x10-6 
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Table 2 

Calculated KIEL for 600 MeV per nucleon ions 

ON CARBON TARGETS 


Incident Ion 

Primary Ion NIEL 
(MeV*cm 2 /g) 

Primary Fragment NIEL 
(MeV»cm 2 /g) 

12 c 

4.83X10- 3 

1.55x10-6 

16 o 

9.54x10-3 

3.42x10-6 

40 Ar 

6.63x1 O' 2 

1 .00x1 O' 5 

56 Fe 

1 .55x1 0" 1 

1.46X10- 5 


The contribution of primary ion fragmentation to NIEL has 
been shown to be negligible. However, this is not necessarily 
true for target fragmentation. When the incident particle is a 
proton or helium ion, for example, it is known that the nuclear 
non-elastic contribution to NIEL is important at high 
energies. In the following, we will show the general trends 
observed as a function of incident ion atomic number and 
energy with and without our estimate of the contribution of 
target fragmentation. 

Fig. 2 shows results for NIEL due to Coulombic 
interactions only, i.e., it excludes target fragmentation. 
Results for heavy ions with Z-values ranging from I to 26 are 
shown for a silicon target. In the space environment the flux 
of ions with Z greater than 26 (iron) falls off rapidly with 
increasing Z. It is seen that for the energies considered, there 
is a monotonic increase in NIEL with incident ion atomic 
number that is proportional to Z2 to first order, as expected. 



Incident Ion Atomic Number 

Fig 2. Heavy ion NIEL in silicon at high energies, excluding target 
fragmentation. 


Fig. 3 shows results analogous to those in Fig. 2 except that 
the effect of target fragmentation is included. Results for this 
nuclear non-elastic contribution were obtained with the 
MCNPX code for protons and alpha particles [4,5], while 
results for heavier ions were scaled according to the reaction 
cross sections given in the work of Tripathi [14,15], as 


discussed above. Fig. 3 thus fills in what has previously been 
a large gap in our knowledge of NIEL, shown in Fig. 1 by the 
relativistic heavy ion region. When target fragmentation is 
included for the energies considered in this work, the effects 
are most pronounced for protons and alpha particles. W’ith 
increasing incident ion atomic number, the target 
fragmentation becomes less significant because the Coulomb 
cross section rises much more rapidly with Z than does the 
reaction cross section. A rule of thumb is that for the 
materials and energies considered here, target fragmentation 
should be accounted for when incident ion Z-values are less 
than about 10. 



Incident Ion Atomic Number 

Fig 3. Heavy ion NIEL in silicon at high energies, including target 
fragmentation 


Figs. 4 and 5 show the full NIEL results for GaAs and C 
targets, respectively, including target fragmentation. These 
targets have significantly different atomic numbers than Si so 
that a range of target Z-values is covered. The results for 
carbon can also be combined with the results for silicon to 
obtain NIEL values for SiC, a material sometimes discussed 
for' possible use in severe environments. 

Heavy Ion NIEL in GaAs at High Energies 



Incident Ion Atomic Number 
Fig 4. Heavy ion NIEL in GaAs at high energies. 
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Heavy Ion NIEL in Carbon at High Energies 



Incident Ion Atomic Number 
Fig 5. Heavy ion NIEL in carbon at high energies. 

The description of heavy ion NIEL for a broad range of 
energies useful for space applications is now reasonably 
complete. Linear Energy Transfer (LET) values, commonly 
used for studies of ionization effects, are known over a similar 
broad range of energies. Since NIEL is the displacement 
damage analog of LET [16], it is interesting to compare the 
two. This is done in Figs 6 and 7 for carbon and iron ions 
incident on silicon. This allows for quantitative comparisons 
to be made about the mechanisms of the ion’s energy loss due 
to displacements and ionization. 



Fig 6. Comparison of NIEL and LET values for carbon ions incident on 
silicon. 



Incident Iron Ion Energy (MeV) 

Fig 7. Comparison of NIEL and LET values for iron ions incident on silicon. 

It is interesting to review some of the features of these new 
NIEL results in the context of previous results. Table 3 
provides a bird’s eye view of NIEL in silicon as a function of 
incident ion energy. The ions considered are H, He, C, Si and 
Fe. NIEL is characterized with respect to the ion’s threshold 
energy (E Ih ), the peak value of NIEL, the ion energy 

corresponding to the peak NIEL (E pk ), the plateau value of 

the non-elastic NIEL at high energies (NIEL^), and the 
estimated ion energy at which the Coulombic and non-elastic 
contributions to NIEL are equal ( E eq ). 


Table 3 

Features of NIEL in Silicon 


Ion 

E if , 

(eV) 

Peak NIEL 
(MeV*cm 2 /g) 

(MeV) 

NIELhp 

(MeV*crrr/g) 

E eq 

(MeV) 

H 

158 

5.032 

0.001 

0.00116 

70 

He 

48 

54.38 

0.001 

0.00171 

900 

C 

25 

385.4 

0.003 

0.00381 

1.5 x 10 4 

Si 

21 

1362 

0.010 

0.0051 1 

2.0 x 10 5 

Fe 

24 

3000 

0.030 

0.00662 

1.0 x 10 6 


The ion’s threshold energy for creating a displacement is 
based on the expression for the maximum energy transfer to a 
silicon atom by the ion, and the assumption that the threshold 
for atomic displacement is 21 eV. The minimum E th occurs 

when the ion is the same mass as the target. 

As the incident ion energies increase, a maximum NIEL 
value is observed as is shown in the third column. Note that 
the peak NIEL for iron is about 600 times that of protons. The 
peak value always occurs in the region governed by the 
screened Coulomb potential [3] as is evident by examining 
the relatively low energies ( E pk ) at the peak. The new results 

in this paper begin at energies of 200 MeV per nucleon and 
are well into the relativistic range. No significant Coulomb 
screening effect is present at these energies. 
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The plateau value of the non-elastic NIEL at high energy is 
shown in the second column from the right. The values for H 
and He are obtained with the MCNPX code [4,5]. Those for 
the heavier ions are estimated using reaction cross section 
values [14,15] to scale the MCNPX results to the high energy 
region. Finally, the last column gives an estimate of the ion 
energy where the Coulombic NIEL is approximately equal to 
the non-elastic NIEL. The reason E eq increases so rapidly 

with ion atomic number is that the Coulombic contribution 
increases rapidly with Z, as can be seen from the peak NIEL 
values. With increasing energy, the Coulombic contribution 
for heavy ions such as Fe takes much longer to come down to 
values comparable to the nuclear non-elastic values. 

Thus, it appears that for ions with Z-values heavier than 
about 10, the contribution of target fragmentation is a small 
correction to the total NIEL value for most energies of 
interest. 

IV. Conclusions 

Estimates of NIEL for high-energy heavy ions have been 
obtained. The relativistic Coulomb interaction of the incident 
ion with the target is the most important contribution. A 
correction for target fragmentation should also be added, 
which is most significant for high-energy, low-Z ions. The 
effect of fragmentation of the incident ion can be ignored. 
The practical consequences of this is that we can now make 
estimates of NIEL values for all the ions present in cosmic 
rays that have significance for radiation effects applications. 

An important implication of the above conclusion is that 
NIEL spectra can now be constructed which include all 
cosmic ray species of interest. These would play a role similar 
to the LET spectra originally devised by Heinrich [17] for 
studying ionization phenomena. For some applications, this 
may considerably simplify the assessment of displacement 
damage produced in the space environment. 

General methods for calculating NIEL for arbitrary ions, 
ion energies and materials have now been developed [3-5, this 
work]. These will be available in a Windows compatible code 
to be distributed by the NASA Space Environments and 
Effects (SEE) Program. 
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